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t  '  r 

thes'e  conditions.  Asymmetry,  .  according  to  this  research,  t  was 
characterized  by  the  change  in 'the  topological  structure  near 
the  saddle  point  associated  with  the  primary  vortices.  In  ad¬ 
dition,  the ''creation  of  vortex;  by  <the  pinching  of  another  vor- 
tex; was  found  to  occur.  This  phenomenon  was  observed  in  both 
~t>he  ''primary  and  secondary  flow.  (■  f 

The  singular  points  obtained  from  the  two-dimensional  par¬ 
ticle  traces  satisfied  the  topological  summation  rule  for  two- 
dimensional  sections  of  a  three-dimensional  flow.  The  primary 
separation  and  associated  vortices  agree  qualitatively  with  ex¬ 
perimental  data.  The  secondary  separation  region,  however, 
contained  additional  topological  structures  that  have  not  been 
resolved  in  previous  experimental  studies. 
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This  paper  was  prepared  by  Kirk  J.  Vanden  of  the  Computational 
Fluid  Dynamics  Section,  Aerodynamics  Branch,  Aeromechanics  Divi¬ 
sion,  Air  Force  Armament  Laboratory,  Eglin  Air  Force  Base,  Flo¬ 
rida.  The  work  was  performed  under  Work  Unit  2307E126  during  the 
period  from  January  1989  to  September  1989. 

This  paper  presents  an  investigation  of  asyimetric  vortical 
flow  on  a  tangent  ogive  body  based  upon  solutions  obtained  by  a 
thin-layer  Navier-Stokes  code  in  use  in  the  Computational  Fluid 
Dynamics  Section. 
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CHAPTER  I 


INTRODUCTICN 

The  purpose  of  this  research  is  to  investigate  asymmetric  vortical 
flew  around  a  body  of  revolution  at  high  angles  of  attack.  The  phenomena 
of  interest  are  the  sheets  of  vorticity  formed  on  the  leeward  side  at 
moderate  angles  of  attack.  At  low  angles  of  attack  the  flow  separates 
synmetrically  on  botl  jjdes  of  a  body  of  revolution  at  sane  point  between 
the  windward  and  leeward  sides.  The  separated  flow  rolls  up  into  a  pair 
of  syrrmetric  vortices.  At  moderate  angles  of  attack  the  vortices  became 
asymmetric  but  steady  and  move  away  fran  the  body  with  respect  to  the 
downstream  direction,  while  new  vortices  form  underneath  them.  At  even 
higher  angles  of  attack  the  flow  in  the  wake  becomes  unsteady  and  the 
vortices  alternately  shed  off  the  body.  Asymmetric  vortices  can  induce 
large  asymmetric  loads  on  the  body,  even  at  zero  sideslip.  The  vortex 

induced  side  force  can  easily  exceed  the  normal  force1 .  Wind  tunnel 

9 

tests  of  the  F-1 1 1  have  shown  the  vortex- induced  yawing  moment  can 
exceed  the  available  control  authority,  through  full  rudder  deflection, 
by  an  order  of  magnitude.  Also,  vehicle  roll  can  cause  the  vortex 
induced  forces  and  moments  to  switch  from  left  to  right  suddenly, 
especially  on  missiles  flying  uncoordinated  manuevers. 

In  this  analysis,  a  modified  version  of  the  Program  EAGLE  -  Flow 

Solver-1  is  used  to  solve  the  thin- layer  approximation  to  the 

Navler -Stokes  equations  to  obtain  the  separated  flow  solutions.  Program 


1 
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EAGLE  -  Grid  Generator  is  used  to  constrict  the  algebraic  grid.  This 
research  uses  a  framework  for  the  observation  of  aerodynamic  flows  which 

stresses  physical  features  and  descriptive  mathematics4,  lb  provide 
a  systematic  and  rigorous  base  to  this  framework  the  mathematics  of 
topology  is  used.  Ibpology  had  its  beginning  in  the  observation  that 
there  are  a  number  of  geometrical  properties  that  are  unchanged  by 
continuous  deformations.  The  use  of  topology  in  aerodynamics  was  first 

introduced  systematically  by  Legendre^,  who  noted  that  the  important 
structures  of  the  flew  can  be  described  by  a  finite  number  of  singular 
(critical)  points  which  must  conform  to  summation  rules.  The  topology  of 
singular  points  has  proved  to  be  a  very  useful  tool  in  analyzing  skin- 
friction  patterns  and  streamlines  frcm  experimental  studies  or  numerical 
calculations . 

Experimental  studies  of  asymmetric  vortex  flew  have  been  conducted 

by  many  researchers.  Lament^  performed  low  Mach  number  wind  tunnel 
tests  of  an  ogive-cylinder  and  identified  two  mechanisms  for  producing 
asymmetric  flow.  One  mechanism  operates  in  the  fully  laminar  or  fully 
turbulent  regime  and  is  qualitatively  described  by  the  impulsive  flow 
analogy.  The  second  mechanism  observed  was  in  transitional  flew  and  is 
analogous  to  the  two-dimensional  crossflow  on  a  circular  cylinder  in  the 

transitional  separation  regime.  Yanta7  et  al.  performed  wind  tunnel 
tests  of  a  tangent  ogive  in  low  turbulence  incompressible  flow  with 
asyrrmetric  vortices.  These  results  supported  hydrodynamic  instability  as 
being  a  cause  of  vortex  asyrrmetry.  The  leeside  flow  was  examined  in 
detail  including  the  streamlines  and  singular  points  on  two-dimensional 
cross  sections  of  the  flow.  The  effect  of  nose  bluntness  was  found  to 
have  no  effect  on  the  basic  qualitative  topological  structure  in  the 
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cross  flew  plane.  Keener8  conducted  subsonic  and  supersonic  wind 
tunnel  tests  of  ogival,  parabolic,  and  conical  forebodies  over  a  wide 
range  of  Reynolds  numbers.  The  principle  cause  of  vortex  asymmetry  was 
also  found  to  be  the  hydrodynamic  instability  of  the  inviscid  flow  field 
with  additional  asymmetry  being  observed  in  the  boundary  layer  at 
transitional  conditions.  Nose  bluntness  was  observed  to  add  additional 
structure  in  one  case.  An  additional  vortex  pair,  not  observed  on  sharp 
nosed  ogives,  was  observed  above  the  primary  vortices  on  a  blunted  ogive 
forebody.  These  vortices  were  believed  to  be  generated  at  the  nose. 

Chapman4  discusses  the  topology  of  bodies  of  revolution  in  general 
and  notes  that,  except  for  crossflow  shocks,  the  topology  is  basically 

the  same  for  both  subsonic  and  supersonic  flows.  Keener9  et  al.  and 

Keener  and  Chapman10  have  conducted  wind  tunnel  measurements 
of  asymmetric  induced  forces  and  moments  on  various  bodies  of  revolution 
for  both  subsonic  and  supersonic  flows,  with  the  former  study  noting  that 
as  the  Mach  number  approaches  supersonic  values  the  vortex-induced  side 
force  approaches  zero.  The  control  of  vortex  orientation  has  been 

investigated  by  Peake11  et  al.  using  injection  of  air  near  the 

nose,  and  by  Fidler1^  who  used  rotation  of  the  nose  region. 

Several  numerical  investigations  have  also  been  done. 

Ying1 ^  et  al.  conducted  transonic  laminar  thin- layer 
Navier -Stokes  calculations  of  the  asymmetric  vortical  flow  on  a 
hemisphere  cylinder  and  studied  the  connection  between  surface  flow 
topology  and  the  surfaces  of  separation.  More  recently,  Siclari  and 

Marconi1^  calculated  supersonic  laminar  formation  of  steady 
asymmetric  vortices  on  a  cone.  These  calculations  were  converged  to 
machine  zero  using  a  Navier -Stokes  algorithm  modified  with  the  assumption 


4 


of  conical  flow. 

In  this  thesis,  details  of  the  secondary  flew  associated  with 
asynrretric  vortex  flow  are  nunerically  investigated  and  found  to  agree 
qualitatively  with  previous  experimental  results.  In  addition,  the 
cxo6sf lew  topology  for  a  nutter  of  angles  of  attack  is  given  to 
illustrate  the  difference  between  syirmstric  flew  and  various  asymmetric 
flows. 

A  description  of  the  governing  equations  and  the  numerical  algorithm 
is  given  in  Chapter  II  followed  by  a  description  of  the  use  of  topology 
in  the  analysis  of  fluid  flows  in  Chapter  III.  The  numerical  grids  and 
convergence  history  are  given  in  Chapter  IV.  Finally,  Chapter  V  presents 
and  discusses  the  results  of  the  calulations  performed. 


CHAPTER  II 


FLOW  SOLVER 

The  boundary- layer  equations  are  useful  for  many  different  types  of 
viscous  flews.  However,  there  are  many  problems  in  which  the  boundary- 
layer  assunption  is  not  valid  and  a  more  exact  form  of  the  Navier-Stokes 
equations  must  be  used.  The  current  problem  of  the  high  angle-of -attack 
flow  over  a  body  of  revolution  is  an  exanple.  The  thin-layer 
approximation  to  the  complete  Navier-Stokes  equations  is  particularly 
useful  in  this  case  since  separated  and  reverse  flow  regions  can  still  be 
exmputed  but  with  a  significant  decrease  in  ccnputational  requirements. 

In  the  formal  derivation  of  the  unsteady  boundary-layer  equations  all 
terms  on  the  order  of  1/(ReL  ),/2  and  smaller  are  neglected.  This 
results  in  the  viscous  terms  containing  derivatives  in  the  directions 
parallel  to  the  body  being  dropped  since  they  are  significantly  smaller 
than  the  viscous  term  normal  to  the  boundary.  In  addition,  the  momentum 
equation  is  reduced  to  a  simpler  equation  in  which  the  pressure  normal  to 
the  wall  is  aonstant.  In  the  thin-layer  approximation  the  viscous  terms 
parallel  to  the  body  are  dropped  again,  but  all  other  terms  of  the 
momentum  equation  are  retained.  The  thin- layer  Navier-Stokes  equations 
are  a  mixed  set  of  hyperbolic-parabolic  partial  differential  equations. 

The  thin-layer  Navier-Stokes  Program  EAGLE  flow  solver  calculates  the 
viscous  terms  explicitly  and  implements  them  after  the  inviscid  terms 
have  been  calculated  in  an  implicit  upwind  method.  The  inviscid  portion 
of  the  algorithm  optionally  uses  either  flux-vector  split  or  flux- 
difference  split  finite-volume  formulations.  By  computing  the  viscous 
terms  explicitly  the  cost  of  the  thin- layer  algorithm  is  only  about 
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15%  greater  than  the  Euler  algorithm  on  an  identical  grid.  This  is  due 
to  the  separate  calculation  of  the  inviscid  terms  for  which  upwind 

differencing  is  most  effective.  Whitfield1  6  developed  the 

baseline  algorithm  with  flux-vector  splitting,  Be  Ik  and  Whitfield16 
extended  it  to  obtain  time-accurate  solutions  on  nulti-blodc  grids,  and 

Gatlin  and  Whitfield17  added  the  thin-layer  approximation. 

The  version  of  the  Program  EAGLE  flow  solver  used  in  this  study  was 
further  modified  by  Sirrp6on18  to  include  a  flux-difference 

splitting  scheme  due  to  Roe1 9.  This  scheme  was  used  for  all 
calculations  performed. 

2.1  Thin-Layer  Navier -Stokes  Approximation 

In  the  thin-layer  approximation  all  viscous  terms  except  those 
normal  to  the  body  are  eliminated  from  the  full  Navier-Stokes  equations. 
The  thin-layer  Navier-Stokes  equations  are  transformed  using  generalized 
boundary  conforming  curvilinear  coordinates  to  allow  the  implementation 
of  boundary  conditions  for  arbitrary  geometries.  These  curvilinear 
coordinates  are 


5  =  £(x,y,z) 

T]  •  ij(x,y,z)  (1 ) 

f  =  £(x,y,z) 

T  =  t. 


Applying  this  transformation,  the  non-dimensional,  divergence  form 
of  the  time-dependent  thin- layer  Navier-Stokes  equations  become 

9Q  +  9£  +  3G  +  3H.a^ 

0T  05  01)  0C  011  ' 


(2) 
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where 


P 

Pu 

Q  =  J  pv 
pv 


F  =  J 


G  =  J 


puU  +  ^xP 
pvO  +  ?yP 

P^  + 

U(etP) 

> 

puV  +  |jxP 

PW  +  Tlyp 
pwV  +  T)2P 
v(efP) 


andU  =  F  u  +  f  v  +  P  w, 
*»x  *»y  *»  z 


and  V  = 


lxU  +  V 


H  =  J 


PW 

puW  +  CxP 
PVW  -■  CyP 
pWW  +  CzP 
w(e*P) 


andw  =  r  u  +  r  v  +  rw, 
» x  » y  >>  z 


T1xtxx  +  *lyTxy  +  *lzTxz 

**  “  HxT yx  +  TlyTyy  +  ^z^yz 

TlxTzx  +  TlyTzy  +  *lzTzz 

(1tx  '“Txx  *  vrxy  +  "Txz  -  Sc1  + 

1y  (uTyx  +  ^yy  +  "Tyz  ’  V  + 

Hz  (uT zx  +  ^zy  +  WIzz  -  qz’  1  _  • 

Here  u,  v,  and  w  are  the  velocity  components  in  the  x,  y,  and  z 
cartesian  coordinates,  respectively,  p  is  the  density  of  the  gas,  P  is 
the  pressure,  and  y  is  the  ratio  of  the  specific  heats.  U,V,  and  W 
are  the  contra var iant  velocities  and  J  is  the  Jacobian  of  the  inverse 


8 


transformation  and  is  given  by 

J '  Vyn\-vc  1 '  VVc'Vc 1  +  Wc^Vc  K  <4) 

Assuming  a  perfect  gas  the  total  specific  energy  of  the  fluid  can  be 
given  as 

e  =  [P/(y-1 )]  +  (p/2)[u2  +  v2  +  v2].  (3) 

Ebr  this  finite  volume  method  the  Jacobian  is  the  volume  of  a 
computational  cell.  The  metric  quantities  are  given  by 


£  ■  J-1  (y  z  -z  y 

Sx  T1  C  ^  C 


V'  fix  ■  ^Vr'Yc’’ r*  ’  r'(y5!,iTiil 


£  =  J_1(z 

Sv 


%  ■  ^'YcYV'  =  ^‘WYt, 


tz  -  J"(Vc‘yn’V1'  =  J",<y51YYr’'  r*  *  ■'"YirYV 


and  give  the  projected  area  of  cell  faces  on  the  cartesian  coordinate 
planes  when  multiplied  by  the  Jacobian.  The  stress  terms  are 


Txx  “  3  He°°  ^2t1xU7}  ~  TlyVT)  ~  ^z*!)^ 

Tyy  3  3  Ke°°  *-2l1yVT)  ”  ^lxUT)  ~  ^z'V 

Tzz  "  3  Re00  *-2TlzWn  ”  ^xUT)  ~  ^y^ 

id 

T*y  "  Re°°  ^yUl)  +  ^x^ 

Txz  Re00  ^l)zUi)  +  MXW1)^ 

ii  M 

T  =  u  —  [n  W  +  n  V  ]  . 

Tzy  Re  L'ly  1)  Mz  ir 
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The  heat  flux  terms  in  the  thin-layer  assumption  are 


2.2  Boundary  Conditions 

If  the  governing  equations  are  cast  in  characteristic  variable  form 
there  will  be  an  eigenvalue  associated  with  each  characteristic  variable 
indicating  the  direction  in  which  information  is  prorogated  across  each 
constant  computational  surface.  This  is  the  basis  for  determining 
boundary  conditions  known  as  characteristic  variable  boundary  conditions. 
This  technique  is  implemented  by  a  layer  of  phantom  points  surrounding 
the  block  boundaries  at  which  the  change  in  dependent  variables  is  taken 
to  be  zero.  Since  these  boundary  conditions  are  treated  explicitly,  they 
are  first-order  accurate.  Characteristic  variable  boundary  conditions 
are  employed  at  all  far  fie  Id  boundaries  except  the  one  at  the  rear  of  the 
body.  It  was  found  that  when  using  characteristic  variable  boundary 
conditions  on  this  boundary,  negative  densities  were  generated  at  the 
phantom  points  and  propagated  forward  through  the  subsonic  boundary 
layer.  Therefore,  the  flow  solver  was  modified  to  implement  extrapolated 
boundary  conditions  in  this  region. 

No-slip  conditions  are  applied  to  the  momentum  terms  at  impermeable 
surfaces.  The  momentum  values  at  the  body  surface  are  mirrored  with 
negative  values  at  the  image  point.  This  gives  a  zero  momentum  flux  at 
solid  wall  boundaries.  The  normal  temperature  derivative  is  set  to  zero 
at  the  shearing  surface  for  adiabatic  wall  conditions.  Since  the  normal 


pressure  gradient  is  assumed  to  be  negligible  at  the  solid  wall,  the 
boundary  conditions  for  energy  and  density  are  given  as 


CHAPTER  III 


TOPOLOGY 

The  mathematical  basis  for  the  behavior  of  the  elementary 
singular  points  of  skin-friction  lines  and  the  topological  rules 

these  points  obey  have  been  known  for  a  long  time.  Davey20  and 

Lightill21  were  the  first  to  point  out  that  skin- friction  lines 
must  obey  a  specific  topological  rule:  the  sum  of  nodal  points 
and  foci  must  exceed  the  number  of  saddle  points  by  two.  Nodal 
points,  foci,  and  saddle  points  are  examples  of  singular  points 
(Figure  1 ).  Singular  points  are  divided  into  two  types  known  as 
nodes  and  saddle  points.  Node  points  are  further  divided  into 
nodal  points  and  foci.  These  singular  points  in  the  skin- 
friction  lines  occur  at  isolated  locations  on  the  surface  where 
the  surface  vorticity  and  skin- friction  become  identically  zero. 

A  nodal  point  is  the  point  that  is  common  to  an  infinite 
number  of  skin-friction  lines.  All  the  lines  associated  with 
that  point,  except  one,  are  tangential  to  a  single  common  line.  At 
a  nodal  point  of  separation,  all  the  skin  friction  lines  are 
directed  inward  toward  the  node.  At  a  nodal  point  of  attachment, 
all  the  skin-friction  lines  are  directed  away  from  the  node.  A 
focus  differs  from  a  nodal  point  in  that  an  infinite  number  of 
skin-friction  lines  spiral  around  the  singular  points.  The  lines 
spiral  away  from  an  attachment  focus  while  a  separation  focus  has 
all  the  lines  spiraling  into  it.  The  singular  point  associated 
with  a  saddle  point  has  only  two  skin-friction  lines  that  pass 
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through  it.  The  halves  of  one  line  are  directed  inwards  towards 
the  center  while  the  halves  of  the  retaining  line  are  directed 
outwards.  These  two  lines  act  as  barriers,  making  one  set  of 
skin-friction  lines  inaccessible  to  an  adjacent  set. 

3.1  Tbpology  of  Streamlines 

J.H.B.  Shtith22,  Perry  and  Fairlie23,  and  J.C.R  Hunt24  et 
al.  were  the  first  to  report  that  the  rules  governing  skin-friction 
line  behavior  are  easily  extended  to  yield  similar  rules 
governing  behavior  of  the  flow  itself.  The  singular  points  of 
streamlines  in  any  of  the  sectional  views  of  a  steady  three- 
dimensional  flow  must  obey  topological  rules  also.  All  of  the 
singular  points  that  are  observed  in  skin-friction  lines  also 
appear  in  the  topology  of  streamlines.  In  addition,  half -nodes 
and  half-saddles  can  occur  on  the  boundary.  A  singular  point  in 
these  sectional  views  is  a  point  (x  ,y  )  where  the  velocity 

is  identically  zero.  This  can  be  stated  as 


u(^,yo)  =  0  , 
v(x^,yQ)  =  0  . 


(6) 


With  the  assurption  of  a  continuous  vector  field,  isolated  critical  points 
conform  to  topological  surmation  rules.  For  a  three-dimensional 
siiiply  connected  body,  the  critical  points  in  a  sectional  view 
must  satisfy  the  relation 


(JN  +  )  -  (  ss  +  )  =  -1 


(7) 
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where  M  denotes  nodes  and  S  indicates  saddles  with  the  primes  denoting 
half  singular  points. 

3.2  Topology  of  Separation 

Topology  has  also  played  a  key  role  in  defining  the  nature 

of  separated  flows.  Tbfcak  and  Peake^  have  given  a 
definition  of  attached  and  separated  flews  in  which  they  define 
two  types  of  separation  called  global  and  local.  An  attached 
flew  is  defined  as  a  flow  containing  two  and  only  two  singular 
points  in  the  skin-friction  lines.  These  two  points  are  a  nodal 
point  of  attachment  and  a  nodal  point  of  separation .  Global 
separation  occurs  when  the  skin-friction  patterns  contain  more 
than  two  nodal  points.  For  every  nodal  point  over  two  a  saddle 
point  must  exist.  This  topological  rule  is  stated  as 

JN  -  JS  =  2  (8) 

For  example,  when  the  flow  contains  three  nodal  points  of  which 
two  must  be  the  same  type  the  lines  emanating  from  these  alike 
nodal  points  must  be  prevented  from  crossing.  The  separatrixes 
of  a  saddle  point  provide  the  needed  barrier.  A  separation  line 
is  defined  as  one  of  these  two  separatrixes.  Tb  summarize  Tbfcak 

and  Peake2 skin-friction  lines  for  attached  flew  cannot 
contain  any  saddle  points,  while  they  must  contain  at  least  one 
saddle  point  in  globally  separated  flow. 

Crossflow  separation,  however,  does  not  always  obey  this 
topological  rule.  The  skin-friction  lines  on  a  body  with  only  a 
cross flew  type  of  separation  contain  only  lines  connecting  the 
nodal  point  of  attachment  at  the  nose  with  the  nodal  point  of 
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separation  at  the  rear.  This  is  the  sane  topological  pattern 
that  exists  without  crossflow  separation .  Even  though  no 
separation  saddle  point  exists,  converging  streamlines  can 
squeeze  together  at  a  rate  that  forces  the  fluid  away  from  the 
body.  The  line  on  which  the  limiting  streamlines  eventually 
converge  appears  to  emerge  from  the  region  of  the  attachment  node 
on  the  nose.  This  type  of  separation  is  known  as  local 
separation.  The  concept  of  limiting  streamlines,  first  used  by 

Sears26  and  later  extended  by  Legendre27,  becomes  important 
for  determining  the  point  at  which  the  flow  is  forced  away  from 
the  body.  Limiting  streamlines  are  streamlines  that,  except  near 
flow  separation,  exist  infinitesimally  close  to  the  surface  of  the 
body  and  their  direction  is  the  same  as  the  skin-friction  lines. 
Therefore,  the  limiting  streamlines  can  be  viewed  as  a  surface  of 
streamlines  whose  projection  on  the  body  are  skin-friction  lines. 

Limiting  streamlines  rising  on  either  side  of  the  line  of 
separation  are  prevented  from  crossing  by  a  stream  surface  which 
emanates  from  the  line  of  separation  itself.  On  a  body  of 
revolution  the  dividing  surfaces  usually  roll  up  into  vortices  as 
they  progress  downstream. 


CHAPTER  IV 


NUMERICAL  CALCULATIONS 

The  configuration  selected  for  this  research  is  a  five 
caliber  tangent  ogive  nose  with  a  circular  afterbody  eight 
diameters  long.  This  configuration  was  selected  because  of  the 
wide  range  of  experimental  studies  performed  using  tangent  ogive 
shapes. 

4.1  Numerical  Grids 

The  grids  used  in  this  analysis  were  generated  using  the 
Program  EAGLE  -  Numerical  Grid  Generation  System28 
developed  jointly  by  the  Air  Force  Armament  Laboratory's 
Ccnputational  Fluid  Dynamics  Section  and  the  Department  of 
Aerospace  Engineering  of  Mississippi  State  University. 

The  single  block  algebraic  fine  grid  (Figures  2  and  3)  has 
70  axial,  60  normal,  and  100  circumferential  points.  This 
equates  to  a  total  of  420,000  grid  points.  The  thin- layer 
Navier-Stokes  solutions  on  this  grid  required  over  1 1 4  million 
words  of  core  memory  on  the  Cray-2  machine  at  the  United  States 
Air  Force  Supercomputer  Center -Kirt  land  ( AESCC-K ) .  This  grid  was 
used  for  all  calculations  except  for  a  grid  refinement  study.  A 
single  block  algebraic  grid  (Figures  4  and  5)  with  60  axial,  50 
normal,  and  61  circumferential  points  was  used  for  a  grid 
refinement  study  at  30.0  degrees  angle  of  attack.  The  results  of 
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Figure  3«  Algebraic  Fine  Grid  -  Back  Plane 
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Figure  4.  Algebraic  Coarse  Grid  -  Side  View 
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Figure  5.  Algebraic  Coarse  Grid  -  Back  Plane 
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the  grid  refinement  study  are  discussed  in  the  analysis  chapter. 
4.2  Convergence  History 

Due  to  the  slow  rate  of  convergence  typical  of  viscous 
solution  algorithms,  it  was  Impractical  to  obtain  solutions 
converged  to  machine  zero.  The  L2  norm  of  the  residual  were 
reduced  between  three  and  four  orders  of  magnitude  for  all 
calculations.  1b  insure  that  the  flew  structure  was  not 
dependent  on  the  convergence  level,  contours  of  flew  variables 
were  oenpared  to  insure  that  no  visible  change  had  oocured  in  the 
last  500  iterations .  Typical  convergence  histories  are  shew.,  in 
Figures  6  and  7.  Figure  6  is  a  symmetric  vortex  case  at  an  angle 
of  attack  of  25  degrees,  and  Figure  7  shews  the  convergence 
history  for  a  30  degree  angle  of  attack  case,  which  is  an 
asynmetric  condition. 


CONVERGENCE  HISTORY 

MACB=1 . 40  RE  =  2 . 0E05  ALPHA=25. 
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Figure  6.  Convergence  History  -  Alpha  z  25.0  deg 


CONVERGENCE  HISTORY 

MACH=1.40  RE=2 -  0E05  ALPHA=30 . 


23 


(W80N21  )001 


0.00  250.00  500.00  750.00  1000.00  1250.00 

ITERATION 

Figure  7.  Convergence  History  -  Alpha  =  30. 0  deg 


CHAPTER  V 


ANALYSIS 

Flow  solutions  are  presented  for  the  five  caliber  tangent  ogive 
configuration  at  supersonic  high  angle  of  attack  conditions.  The  Mach 
nunber  was  1 .4  for  all,  calculations  and  the  angle  of  attack  varied  frcm 
25  to  40  degrees.  The  Reynolds  number  was  200,000  based  on  body  base 
diameter.  These  flight  conditions  were  chosen  because  of  the  wide  range 
of  asyrrmetric  vortical  flow  patterns  that  arise.  All  calculations  were 
performed  assuming  laminar  flow  so  as  to  eliminate  any  effects  a 
turbulence  model  may  have  on  the  onset  and  development  of  asynmetry. 

Based  on  Binion  and  Stanewsky29,  the  assumption  of  laminar 

flew  at  this  Reynolds  number  and  high  angles  of  attack  is  a  valid  one. 

In  this  analysis  all  streamlines  were  computed  using  FLOT3D.30 
This  program  uses  a  2nd  order  Runge-Kutta  integration  routine  with  five 
integration  steps  per  grid  cell.  The  streamlines  on  all  of  the  •  > 

two-dimensional  cross  -sect  ions  of  the  flow  were  calculated  on  a  plane  one 
half  cell  f award  of  the  last  constant  i  grid  plane  since  it  is 
perpendicular  to  the  body  longitudinal  direction.  Since  the  plane  of 
cell  centers  is  not  exactly  perpendicular,  the  longitudinal  momentum 
conponent  was  set  to  zero  for  plotting  purposes.  A  non-zero  momentum 
term,  while  using  PLOT3D,  will  give  incorrect  two-dimensional  streamlines 
if  the  grid  is  not  perpendicular  to  the  body.  These  streamline  plots  are 
viewed  from  the  upstream  direction.  Since  these  streamline  calculations 

involve  only  two  of  the  three  velocity  components  they  are  not 
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necessarily  projections  of  the  three-dimensional  streamlines  onto  the 
cross-sectional  plane.  A  particle  injected  into  the  flow  at  a  specific 
location  will  not  neocessarily  follow  the  streamlines  given  by  the 
two-dimensional  calculations.  During  the  analysis,  sufficient 
three-dimensional  streamline  calculations  were  performed  to  Insure  that 
the  two-dimensional  calculations  represent  true  three-dimensional 
structures.  For  example,  three-dimensional  streamlines  originating  near 
the  nose  were  traced  to  the  end  of  the  body  to  insure  that  the  vortex 
carters  are  located  in  the  spots  indicated  by  the  two-dimensional 
calculations  (Figure  8).  Also,  the  points  of  attachment  and  separation  as 
given  by  the  limiting  streamlines  were  carpared  to  the  points  of 
attachment  and  separation  suggested  by  the  two-dimensional  streamlines 
(Figure  9). 

Figure  10  shows  the  streamlines  at  25  degress.  The  flow  is  symmetric 
at  this  angle  of  attack.  This  two-dimensional  streamline  calculation  is 
on  a  cross-section  one  half  cell  forward  of  the  back  of  the  body.  This 
location  is  X/L=97.0%.  The  flow  is  dominated  by  two  large,  symmetric 
vortices  formed  by  the  separation  and  roll  up  of  the  shear  layer.  Figures 
11  and  12  shew  the  left  and  right  secondary  flow  regions,  respectively. 

Each  secondary  flow  region  contains  three  nodes,  1  saddle,  and  four  half 
saddle  points.  Figure  13  shows  the  crossflow  topology  for  this  case. 

This  is  in  agreement  with  equation  (7),  and  thus  is  a  kinematically  possible 
flew. 

Figure  14  shows  the  vortical  flow  pattern  at  30  degrees  angle  of 
attack.  In  this  case  the  flow  is  dominated  by  two  large,  slightly 
asymmetric  vortices.  These  primary  vortices  induce  a  complicated 
structure  of  secondary  separations,  attachments,  and  vortices.  Figure  15 
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Figure  8.  Comparison  of  Three-Dimensional  and  Two-Dimensional  Streamlines 


Figure  9.  Comparison  of  Limiting  Streamlines  and  Two-Dimensional  Streamlines 
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Figure  11.  Left  Secondary  Flow  Streamlines 
Alpha  =  25.0  deg,  X/L  =  97. 0% 
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Figure  12.  Right  Secondary  Flow  Streamlines 
Alpha  =  25.0  deg,  X/L  =  97.0% 
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Figure  13.  Crossflow  Topology  - 
Alpha  =  25.0  deg,  X/L  =  97. 0* 
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shews  the  limiting  streamlines  plotted  on  the  back  portion  of  the  body. 
The  view  shown  is  one  in  which  the  observer  is  at  the  rear  of  the  body 
locking  down  inside  the  body  towards  the  no6e.  The  first  layer  of  grid 
points  is  shewn  to  give  a  sense  of  perspective.  The  blank  region  at  the 
stagnation  point  is  the  block  to  block  connection.  There  are  12  points  of 
separation  and  attachment  visible.  The  primary  separation  node6  at  the 
rear  of  the  body  look  as  if  they  each  consist  of  two  unconnected  groups 

of  streamlines.  These  are  connected  near  the  ncee  as  shown  in  Figure  16 

» 

but  split  apart  about  a  third  of  the  way  down  the  body.  They  briefly 
reconnect  before  splitting  a  second  and  final  time  (Figure  17).  Several 
places  were  observed  where  streamlines  connect  these  two  lines  but  these 
were  very  few.  It  is  not  known  whether  this  is  a  feature  of  the  flow, 
the  effect  of  decreasing  grid  resolution  near  the  rear  of  the  body,  or 
the  lack  of  a  sufficient  concentration  of  particle  traces  generated  by 
the  plotting  package.  Grid  resolution  appears  to  be  the  least  likely 
since  the  two  regions  reconnect  in  a  coarser  region  of  the  grid  than 
where  they  first  separated.  A  two-dimensional  streamline  calculation  in 
a  cross-sectional  view  of  this  region  shows  no  unusual  structure 
associated  with  the  primary  separation  line  (Figure  18).  All  of  the 
separation  and  attachment  lines  start  at  the  nose  (Figure  19). 

A  detailed  streamline  calculation  of  the  left  secondary  separation 
region  is  shown  in  Figure  20.  There  are  three  nodes,  one  saddle,  and 
four  half  saddle  points  in  this  region.  The  two  vortices  closest  to  the 
body  ( Figure  2 1 )  are  associated  with  separation  and  reattachment  half 
saddle  points,  while  the  remaining  vortex  is  associated  with  a  full 
saddle  point  and  is  separated  from  the  other  two  vortices  by  streamlines 
frem  the  primary  vortices.  This  vortex  does  not  appear  to  be  directly 


Figure  16.  Primary  Separation  Limiting  Streamlines 
Nose  Region,  Alpha  =  30. 0  deg 
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Figure  17  Primary  Separation  Limiting  Streamlines 
Mid-Body  Region,  Alpha  =  30. 0  deg 
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Figure  20.  Left  Side  Secondary  Flow  Streamlines 
Alpha  =  30.0  deg,  X/L  =  97.0$ 


Figure  21.  Left  Side  Secondary  Flow  Streamlines 
Alpha  =  30.0  deg,  X/L  =  97. 0* 
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related  to  the  secondary  separations  and  attachments.  Instead,  It 
appears  to  be  formed  when  the  streamlines  from  the  primary  vortex  are 
turned  from  moving  towards  the  windward  side  to  the  leeward  side  by  the 
free  stream  flow.  This  secondary  flow  pattern  was  observed  in 
experimental  two-dimensional  cylinder  flow  by  Bouard  and 

Goutanceau31 .  The  secondary  flow  on  the  other  side  of  the 
cross-section  is  similar  (Figures  22-24). 

There  are,  in  total,  8  foci,  3  saddle  points,  and  12  half  saddle 
points  in  the  entire  cross-sectional  plane.  In  addition  to  the  16 
singular  points  associated  with  the  secondary  separation  regions,  there 
are  a  half  saddle  at  the  windward  stagnation  point,  two  half  saddles 
corresponding  to  the  primary  separations,  two  foci  denoting  the  primary 
vortices,  a  saddle  at  the  top  of  the  primary  vortices,  and  a  half  saddle 
where  the  primary  vortices  induce  attachment  at  the  top  of  the  body 
(Figure  25).  Equation  (7)  can  now  be  used  to  test  the  validity  of  this 
topology.  Substituting  the  values  of  nodal,  half  nodal,  saddle,  and  half 
saddle  points  gives 

(8+^)-(3+4|)  =  -1 

which  is  in  oonplete  agreement.  Therfore,  this  is  a  kinematically 
possible  flow  structure.  The  primary  structure  agrees  with  that  given  by 

Yarrta^  et  al.  for  a  tangent  ogive  configuration  in  irccrpressible 
flew. 

The  grid  refinement  study  was  performed  at  this  angle  of 
attack.  Figure  26  shows  the  two-dimensional  streamlines  for  the  coarse 
grid  at  the  same  location  as  Figure  14  for  the  fine  grid.  The  basic 
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Figure  22  •  Right  side  Secondary  Flow  Streamlines 
Alpha  =  30.0  deg,  X/L  =  97. OX 
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Figure  23*  Right  Side  Secondary  Flow  Streamlines 
Alpha  =  30.0  deg,  X/L  =  97.0$ 


Figure  24.  Right  Side  Secondary  Flow  Streamlines 
Alpha  =  30.0  deg,  X/L  =  97. 0* 
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Figure  25.  Crossflow  Topology  - 
Alpha  =  30.0  deg,  X/L  =  97.05 


46 


47 


structure  of  the  primary  vortices  remained  the  same.  The  secondary  flow 
is  not  the  same,  but  is  similar.  Much  of  the  detail  shown  in  Figures 
20-24  is  missing.  This  is  to  be  expected  since  the  grid  resolution  was 
much  lower  in  this  region.  This  may  indicate  that  if  the  dimensions  of 
the  fine  grid  used  was  increased  substantially,  additional  singular 
points  might  appear.  The  limiting  streamlines  are  shewn  in  Figure  27. 

Flexure  28  shows  the  streamlines,  at  X/L»97.0%,  for  an  angle  of 
attack  of  37.5  degrees.  Calculations  were  not  done  to  see  how  the 
structure  changed  between  30  degrees  (Figure  14)  and  37.5  degrees.  With 
such  a  large  gap  in  angle  of  attack  it  is  possible  that  the  flow 
structure  passed  through  sane  intermediate  step  rather  than  the  change 
portrayed  by  just  looking  at  the  two  cases  7.5  degrees  apart.  Comparing 
30  degrees  (Figure  14)  and  37.5  degrees  (Figure  28),  one  primary  vortex 
has  moved  above  the  other  one  and  is  very  irregular  in  shape.  The  center 
of  this  vortex  is  located  approximately  two  body  diameters  away  from  the 
body.  Figure  29  shows  the  three-dimensional  streamlines  that  form  this 
vortex.  The  other  vortex  is  much  closer  to  the  body  than  in  the  30 
degree  case.  Figure  30  shows  an  enlargement  of  the  left  secondary  flow 
region.  A  new  node  and  saddle  oanbination  has  appeared  within  one  of  the 
previous  secondary  flow  vortices.  This  new,  and  very  small  vortex, 
appeared  when  a  portion  of  the  vortex  containing  it  flattened  and  pinched 
off  (Figure  31 ).  The  right  side  (Figure  32)  retains  the  same  stucture  as 
the  30  degree  case  (Figure  22).  Except  for  the  left  secondary  flow 
region  and  asyrrmetric  primary  vortices,  the  37.5  case  is  very  similar  to 
the  30  degree  case  in  terms  of  location  and  types  of  singular  points. 

With  the  addition  of  a  new  foci  and  saddle,  ccnpared  to  the  30  degrees 
(Figure  25),  the  simulation  equation  is  still  satisfied.  Figure  33  shows 


Figure  28.  Streamlines  -  Alpha  =  37.5  deg,  X/L  =  97  .Of 


Figure  29.  Three-Dimensional  Streamlines  of  Left  Vortex 
Alpha  =  37.5  deg,  X/L  =  97.0$ 
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Figure  31.  Left  Side  Secondary  Flow  Streamlines 
Alpha  s  37-5  deg,  X/L  =  97.0? 
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Figure  32.  Right  Side  Secondary  Flow  Streamlines 
Alpha  =  37.5  deg,  X/L  =  97. Of 


Figure  33.  Rear  Body  Limiting  Streamlines 
Alpha  =  37.5  deg 
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the  limiting  streamlines  on  the  back  portion  of  the  body  for  the  37.5 
degree  case.  The  overall  pattern  is  the  same  as  the  30  degree  case 
(Figure  15)  although  sane  skewness  can  be  seen.  The  attachment  node 
resulting  fran  the  primary  vortices  is  offset  to  the  left  of  the  body 
center  line.  The  other  nodes  are  slightly  asymmetric  as  well. 

Figures  34-36  show  the  density  contours  at  three  axial  locations 
along  the  body  for  37.5  degrees.  Figure  34  is  the  same  plane  as  the 
streamlines  in  Figure  28.  The  stong  gradient  in  density  between  the 
vortices  is  a  crossflow  shock.  The  density  contours  show  what  appears  to 
be  two  circular  structures  where  the  primary  vortices  are  usually  located 
in  crossflow  separation .  However,  this  is  not  the  correct  position  as 
can  be  seen  frcm  oarparison  with  the  streamlines.  The  density  contours 
indicate  that  the  left  vortex  is  closer  to  the  body  than  it  actually  is. 
Figures  37-39  are  constant  vorticity  contours  at  the  same  axial  location 
of  the  density  contours,  again  at  37.5  degrees.  They  also  do  not  give 
the  correct  position  of  the  left  vortex.  These  calculations  indicate 
that  it  is  possible  to  have  nearly  symmetric  vorticity  and  density 
countour  plots  when  there  are  large  asymmetries  in  the  flow. 

A  calculation  was  ocnpleted  at  37.5  degrees  with  the  normal  viscous 
term  set  to  zero  frcm  20  grid  points  away  from  the  body  to  the  far  fie  Id. 
Figure  40  shows  the  streamlines  and  can  be  ocrpared  with  Figure  28. 

Since  the  streamlines  were  basically  unchanged,  the  flow  outside  the 
boundary  layer  region  is  essentially  inviscid  and  viscosity  does  not  play 
a  critical  role  in  the  relative  size,  shape,  and  location  of  the  primary 
vortices.  Based  on  the  small  effect  of  this  normal  viscous  term,  the 
thin-layer  approximation  produced  accurate  results  away  frcm  the 
boundary  layer  region. 


Figure  3H.  Constant  Density  Contours 
Alpha  :  37.5  deg,  X/L  =  97.0* 
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Figure  35.  Constant  Density  Contours 
Alpha  =  37.5  deg,  X/L  =  31. OS 
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Figure  36.  Constant  Density  Contours 
Alpha  =  37.5  deg,  X/l  =  15. 0% 
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Figure  37.  Constant  Vorticity  Contours 
Alpha  =  37.5  deg,  X/L  =  97. 0* 
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Figure  38.  Constant  Vorticity  Contours 
Alpha  =  37.5  deg,  X/L  =  31. 0% 


Figure  39.  Constant  Vorticity  Contours 
Alpha  =  37.5  deg,  X/L  =  15. 0* 
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Fic^ire  41  shows  the  two-dimensional  streamlines  for  38.5  degrees. 

The  additional  degree  from  37.5  degrees  to  38.5  degrees  has  resulted  in  a 
third  primary  vortex  appearing  and  an  additional  saddle  point  appearing 
between  one  of  the  previous  primary  vortices  and  the  new  one.  Figure  42 
shews  an  enlargement  of  the  area  containing  the  new  vortex  and  the  new 
primary  saddle  point.  Figure  43  shews  the  three-dimensional  streamlines 
that  form  the  core  of  the  new  \ortex.  The  primary  vortex  that  was 
located  the  farthest  from  the  body  in  the  37.5  degree  case  has  moved  even 
farther  away  in  the  38.5  degree  case.  These  two  vortices  rotate  in  the 
same  direction  and  are  fed  by  the  same  shear  layer.  They  are  separated 
by  a  separatrix  of  the  saddle  point.  The  vortlcity  contours  in  Figures 
44-49  show  this  shear  layer  at  different  longitudinal  locations.  This 
third  vortex  appears  to  occur  when  the  primary  vortex  flattens 
sufficiently  to  pinch  off  a  portion  of  itself,  thus  forming  a  separate 
region  of  rotation.  This  agrees  well  with  the  conceptual  three  vortex 

configuration  given  by  Peake  and  Tbbak32.  Figures  50  and  51 
shew  the  entire  left  side  secondary  flew  region,  and  Figure  52  shews  the 
right  side.  The  secondary  flows  on  both  sides  contain  the  same  singular 
points  as  the  37.5  degree  case  (Figures  30-32).  The  additional  foci  and 
saddle  point  in  the  primary  flow  region  balance  each  other  in  the 
topological  summation  equation .  In  the  entire  cross-flcw  plane  there  are 
12  half  saddles,  10  foci,  and  5  saddles  (Figure  53).  Substituting  into 
the  topological  summation  equations  gives 

(  10  +  -  (  5  +  4|)  =  -i 

in  ccrrplete  agreement.  The  limiting  streamlines  on  the  rear  of  the  body 


Figure  41.  Streamlines  -  Alpha  =  38.5  deg,  X/L  =  97 
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Figure  42.  Streamlines  -  Alpha  =  38.5  deg,  X/L  =  97 .OX 
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Figure  1|3*  Three-Dimensional  Streamlines  of  New  Vortex 
Alpha  =  38.5  deg,  x/L  =  97. OS 


Figure  44.  Constant  Vorticity  Contours  - 
Alpha  =  38.5  deg,  X/L  =  55. 0* 
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Figure  il5.  Constant  Vorticity  Contours 
Alpha  =  38.5  deg,  X/L  =  52.05 
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Figure  '(6.  Constant  Vorticity  Contours 
Alpha  =  36.5  deg,  X/L  =  50. OX 
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Figure  n 9.  Constant  Vorticity  Contours  - 
Upha  =  38.5  deg,  X/L  =  30.0? 
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Figure  50.  Left  Side  Secondary  Flow  Streamlines 
Alpha  =  ?8.5  deg,  x/L  =  97. 0* 


Figure  51 •  Left  Side  Secondary  Flow  Streamlines 
Alpha  =  38.5  deg,  X/L  =  97. Of 


75 


Figure  52.  Right  Side  Secondary  Flow  Streamlines 
Alpha  =  38.5  deg,  X/L  =  97.0* 
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Figure  53.  Crossflow  Topology  - 
Alpha  =  38.5  deg,  X/L  =  97. 0* 
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are  shown  in  Figure  54.  The  separation  and  attachment  convergence  lines 
are  noticeably  more  skewed  than  the  37.5  degree  case  (Figure  33).  The 
particle  traces  forming  the  first  attachment  line  above  the  primary 
separation  line  on  the  right  side  appear  to  spread  apart.  This  is  from 
the  plotting  package  and  not  part  of  the  flow  structure.  This  software 
draws  particle  traces  in  both  the  upstream  and  downstream  directions. 
Since  the  body  is  terminated,  and  not  infinite,  there  is  a  region  near 
the  base  that  only  particle  traces  from  the  upstream  direction  reach. 
Thus,  the  convergence  lines  formed  by  limiting  streamlines  are  not 
resolved  well. 

Figure  55  shows  the  streamlines  for  40.0  degrees  angle  of  attack. 
Since  no  calculations  were  made  between  38.5  degrees  (Figure  42)  and  40.0 
degrees  it  is  impossible  to  determine  how  the  flow  reached  this  new 
state.  However,  it  is  clear  that  major  changes  took  place  going  from  the 
three  vortex  configuration  of  38.5  degrees  (Figure  42)  to  a  two  vortex 
configuration.  One  vortex  of  the  40.0  degree  calculation  is  very  near 
the  body  and  similar  in  shape  to  the  vortices  observed  at  lower  angles  of 
attack.  The  other  vortex  is  much  larger  and  very  irregular  in  shape. 

The  saddle  point  associated  with  these  two  primary  vortices  is  to  the 
left  of  the  center  of  the  smaller  vortex  and  nearly  in  line  with  the  left 
side  of  the  body.  Figure  56  shows  the  streamlines  in  the  right  secondary 
flew  region.  The  arrangement  and  type  of  singular  points  is  similar  to 
that  already  discussed  for  lower  angles  of  attack  except  a  focl-saddle 
combination  has  disappeared.  The  particle  traces  of  the  left  secondary 
flow  (Figure  57)  show  the  most  complex  topology  yet  calculated.  There 
are  three  foci  and  two  saddle  points  in  addition  to  two  half  saddle 
points  on  the  body.  An  enlargement  of  the  region  with  the  two 
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Figure  54.  Rear  Body  Limiting  Streamlines 
Alpha  =  38.5  deg,  X/L  =  97. Of 
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Figure  57.  Left  Side  Secondary  Flow  Streamlines 
Alpha  =  40.0  deg,  X/L  =  97.0$ 
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foci-sactile  point  combination  is  shewn  in  Figure  58.  This  gives,  in  the 
overall  flow,  eight  nodes,  three  saddles,  and  12  half  saddle  points. 
Substituting  into  the  topology  summation  equation  gives 

(  7+^)-(  2  +  -  -1 

again  in  ocnplete  agreement.  Figure  59  shows  the  limiting  streamlines  on 
the  rear  of  the  body.  They  are  highly  asynmetric  and  contain  a  number  of 
sharp  kinks. 

Figures  60-64  show  the  circumferential  coefficient  of 
pressure  at  X/L=97.0%  for  all  the  angles  of  attack  calculated.  Theta 
equal  to  zero  and  360  correspond  to  the  windward  stagnation  point.  The 
direction  of  rotation  is  clockwise  in  relation  to  the  two-dimensional 
streamline  calculations  previously  shown.  The  symnetric  pressure 

coefficient  distribution  agrees  qualitatively  with  Lament6  and  Ericsson 

and  Reding-*3.  The  buirps  in  the  pressure  coefficient  at 
approximately  a  theta  of  135  degrees  for  each  of  the  cases  correspond  to 
the  portion  of  the  secondary  flow  with  the  small  vortices  nearest  the 
body.  T&ble  1  gives  the  forces  and  moments  for  all  cases  calculated. 
Since  the  free  stream  flow  is  supersonic  the  vortex-induced  side  forces 
sure  smaller  than  the  normal  force,  agreeing  qualitatively  with 

Lament6,  Ericsson  and  Redding34,  and  Keener3  et  al. 
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Figure  58.  Left  Side  Secondary  Flow  Streamlines 
Alpha  =  HO.O  deg,  X/L  =  97.0* 
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Figure  59.  Rear  Body  Limiting  Streamlines 
Alpha  =  40.0  deg,  x/L  =  97.0* 
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Figure  61.  Circumferential  Pressure  Coefficient  - 
Alpha  =  30.0  deg,  X/L  =  97. 0* 


Table  1.  Forces  and  Moments 


CHAPTER  VI 


CONCLUSIONS 


Asymmetric  vortical  flow  has  been  calculated  numerically  using  a 
thin-layer  Navier -Stokes  algorithm.  It  was  found  that  the  thin-layer 

approximation  models  the  physics  acceptably  In  both  the  boundary-layer 

* 

and  outer  flow  regions. 

Both  sytrtnetric  and  asynmetric  cases  were  calculated  and  analyzed. 

The  secondary  flow  was  shown  in  detail.  The  basic  structure  was  found  to 
agree  qualitatively  with  experimental  results,  and  a  number  of  additional 
structures  and  characteristics  were  observed.  The  complexity  of  the 
secondary  flow  region  increased  with  angle  of  attack. 

Crossflow  topology  diagrams  were  given  for  symmetric  and  asymmetric 
conditions  to  illustrate  the  fundamental  difference  between  these 
conditions.  Asymmetry,  according  to  this  research,  was  characterized  by 
the  change  in  the  topological  structure  near  the  saddle  point  associated 
with  the  primary  vortices.  In  addition,  the  creation  of  a  vortex  by  the 
pinching  of  another  vortex  was  found  to  occur.  A  saddle  point  was  formed 
In  the  middle  of  this  two  vortex  structure.  This  phenomenon  was  observed 
in  both  the  primary  and  secondary  flow. 

In  general,  these  results  were  in  qualitative  agreement  with 
experimental  data,  but  showed  additional  structures  not  observed 
previously. 
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